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INTRODUCTION 

Multiple  physical  phenomena  can  be  simulated  with  closed  loop  controls.  There  are 
applications  throughout  the  industry  from  robotics  to  flight  control  systems.  This  can  include  the 
physics  of  thermal,  structural,  and  fluid  flow,  which  can  be  combined  with  software  simulations  of  a 
control  system. 

A  structural  finite  element  analysis  (FEA)  has  been  coupled  with  a  computational  fluid 
dynamic  (CFD)  analysis.  The  FEA  and  CFD  analyses  were  run  simultaneously  using  co-simulation. 
Both  of  these  analyses  interact  with  a  user  subroutine  that  simulates  a  control  system.  In  this 
demonstration  model,  the  canard  is  rotated  (fig.  1)  using  a  piezoelectric  beam  and,  as  the  beam 
bends,  the  canard  rotates.  Voltage  differences  control  the  amount  of  bending  in  the  beam.  The 
voltage  is  controlled  by  a  user  subroutine  that  reads  the  angle  of  the  canard  and  then  varies  the 
voltage  to  reach  the  target  angle.  The  CFD  analysis  applies  a  load  on  the  canard  due  to  the  flow  of 
fluid  around  it. 


Figure  1 

Typical  canards 


MODEL  DESCRIPTION 

The  model  has  three  major  components  (fig.  2),  which  are  a  CFD  portion  to  simulate  fluid 
flow,  a  finite  element  method  portion  to  simulate  structural  response,  and  a  user  subroutine  to 
simulate  a  control  system. 
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Figure  2 

Major  components 


Fluid  Model 

The  CFD  model  consists  of  a  2-D  fluid  domain  with  an  inlet  flow  on  the  left  and  an  exit  on  the 
right  (fig.  3).  There  is  a  boundary  layer  around  the  canard.  The  flow  is  transient,  incompressible,  and 
laminar. 


Flow  Direction 
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Figure  3 
CFD  model 
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The  fluid  structure  interaction  occurs  at  the  interface  between  the  CFD  mesh  and  the  finite 
FEA  mesh  (fig.  4).  There  is  a  two-way  interaction  at  the  interface  that  is  created  by  a  tied  constraint. 
As  the  canard  rotates,  the  mesh  adapts  to  conform  to  the  new  shape. 


CFD  Mesh 


Figure  4 

CFD  and  FEA  interface 


Structural  Model 

The  FEA  structural  model  (fig.  5)  consists  of  a  canard  and  a  piezoelectric  beam.  The  beam 
and  canard  are  connected  by  a  rigid  connector  and  are  both  deformable  bodies.  An  implicit  dynamic 
solver  was  used  to  calculate  the  deformation  of  the  structural  model. 


Figure  5 
FEA  model 
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When  there  is  a  difference  in  voltage  across  the  piezoelectric  material,  the  material  expands 
or  contracts  (fig.  6)  depending  on  the  material  orientation.  This  expansion  and  contraction  causes 
the  beam  to  bend.  As  the  difference  in  voltage  increases,  the  bending  of  the  beam  increases. 


Zero  Voltage 


This  side  expands 


Voltage  Load 


This  side  contracts 


Figure  6 

Operation  of  piezoelectric  beam 


Control  Subroutine 

The  voltage  applied  to  the  piezoelectric  beam  is  controlled  by  a  user  subroutine.  This 
subroutine  can  read  the  current  angle  of  the  canard  and  then  calculate  a  new  voltage  to  be  applied 
to  the  piezoelectric  material.  This  cycle  continues  throughout  the  analysis. 

Two  simple  control  schemes  were  implemented.  In  the  first  case,  the  voltage  is  ramped  up 
quickly  and  then  kept  constant,  and  the  canard  and  beam  are  allowed  to  move  freely.  In  the  second 
case,  the  voltage  was  turned  off  when  the  angle  of  the  canard  exceeded  a  target  value,  and  then  the 
voltage  is  turned  back  on  when  the  angle  of  the  canard  is  below  the  target  value.  These  control 
routines  are  intentionally  simplistic  to  demonstrate  the  coupled  approach.  More  sophisticated 
routines  can  be  implemented  as  needed  by  the  user  by  rewriting  the  FORTRAN  subroutine. 


RESULTS 


Computational  Fluid  Dynamic  Results 

The  pressure  profile  when  the  canard  is  in  its  initial  horizontal  state  is  shown  in  figure  7.  The 
pressure  is  higher  in  the  forward  edge  of  the  canard  and  lower  at  the  trailing  edge  as  expected.  The 
pressure  profile  after  the  canard  has  rotated  is  shown  in  figure  8.  As  expected,  the  low  pressure  area 
has  moved  to  the  top  surface  of  the  canard  and  the  high  pressure  area  has  moved  to  the  lower 
surface  of  the  canard. 
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Figure  7 

Pressure  profile  of  horizontal  canard 
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Figure  8 

Pressure  profile  of  canard  at  an  angle 
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Finite  Element  Analysis  Results 

The  von  Mises  stress  is  shown  in  figures  9  and  10.  In  figure  9,  the  beam  and  canard 
assembly  are  shown  in  their  initial  neutral  state  with  zero  stress.  In  figure  10,  the  assembly  is  shown 
in  its  deformed  state,  where  the  beam  is  bent  and  the  canard  has  rotated.  The  high  stress  in  the 
piezoelectric  material  is  easily  visible  on  the  beam. 
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Figure  9 

Stress  in  piezoelectric  beam  with  zero  voltage 
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Figure  10 

Stress  in  piezoelectric  beam  with  voltage  applied 
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Control  Results 

Plots  of  the  voltage  and  canard  angle  are  shown  in  figures  11,  12,  and  13.  Figure  11  shows 
case  one,  where  the  voltage  is  quickly  ramped  up  and  kept  constant.  The  canard  angle  shows  how 
the  canard  motion  damps  out  overtime  converging  to  a  constant  angle.  Figure  12  shows  case  two, 
where  the  voltage  is  turned  off  when  the  canard  angle  goes  above  the  target  value.  The  voltage  is 
turned  on  again  when  the  angle  drops  below  the  target  value.  The  plot  of  the  canard  angle  shows  it 
is  slowly  damping  down  toward  the  target  angle.  Figure  13  shows  a  comparison  between  the  canard 
rotations  in  case  one  with  constant  voltage  and  case  two  with  voltage  being  turned  on  and  off. 


Figure  1 1 

Plot  of  constant  voltage  and  canard  angle 


Figure  12 

Plot  of  varying  voltage  and  canard  angle 
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Figure  13 

Plot  comparing  canard  angles 


CONCLUSIONS 

A  combination  of  structural  finite  element  analysis  with  computational  fluid  (CFD)  analysis 
controlled  by  a  FORTRAN  user  subroutine  was  successfully  demonstrated.  There  was  a  bidirectional 
interaction  between  all  three  processes  or,  in  other  words,  a  fully  coupled  analysis.  The  CFD  model 
can  be  expanded  to  a  full  3-D  and  third  party  codes,  such  as  STAR-CCM+,  can  be  used  instead  of 
the  ABAQUS  solver.  Both  the  ABAQUS  implicit  and/or  explicit  solver  can  be  used  with  co-simulation. 
The  control  routine  can  be  written  in  both  FORTRAN  and  C++  programming  languages.  This  creates 
a  very  flexible  system  for  solving  complex  multiphysics  problems  with  control  systems. 
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